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Abstract 

Heavy meson pairs produced in the decays of heavy quarkonium resonances at 
e + e~ machines (beauty and tau-charm factories) have the useful property that the 
two mesons are in the CP-correlated states. By tagging one of the mesons as a CP 
eigenstate, a lifetime difference of heavy neutral meson mass eigenstates Ar may 
be determined by measuring the leptonic branching ratio of the other meson. We 
discuss the use of this and related methods both in the case where time dependent 
mixing is small and when it is significant. We consider the impact of possible 
CP-violating effects and present the complete results for CP-entangled decay rates 
with CP-violation taken into account. 
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1 Introduction 



One of the most important motivations for studies of heavy meson mixing is the pos- 
sibility of observing a signal from new physics which can be separated from the one 
generated by the Standard Model (SM) interactions. If H° is a neutral heavy meson, 
AQ = 2 transitions, occurring only at one loop in the Standard Model, as well as possible 
new interactions, generate contributions to the effective operators that change H° state 
into H° state leading to the mass eigenstates 

\H l2 )=p\H°)±q\H ), (1) 

where the complex parameters p and q are obtained from diagonalizing the H° — H° mass 
matrix HJ. It is conventional to parameterize the mass and width splittings between these 
eigenstates by 



m 2 - mi r 2 - Ti 

y = -2f-' (2) 

where and r 12 are the masses and widths of Hi 2 and the mean width and mass 
are Y = (Ti + T 2 )/2 and m = (mi + m 2 )/2. Since y is constructed from the decays of 
H into physical states, it should be dominated by the Standard Model contributions, 
unless new physics significantly modifies AQ = 1 interactions. 

Our goal in this paper is to discuss time independent methods to determine y for 
the various oscillating hadron species (i.e. H = D°, B° and B s ) at electron-positron 
colliders where an HH° pair can be produced in a correlated initial state. In particular 
we will consider methods which are time independent and therefore may be applied in 
experiments where the time history is difficult to obtain. Results which are linear, rather 
then quadratic, in y are of particular interest because y is generically small for B® and 
D° mesons. 

Indeed, in the case of D°, such studies may be carried out at tau-charm factories 
running on the ^(3370) resonance. In this case x is small and it could well be that 
y 3> x. It has been argued that the Standard Model y > x @ with y ~ 0(1%). Thus, 
lifetime differences may be the dominant form of mixing in D° and their study is perhaps 
within the reach of proposed experiments. 

In the case of B%, x = 0.730 ± 0.029 § and y is small, of the order of 0.3% §. 
It is possible that at future B-factories there would be enough statistics to measure y. 
In addition, if e + e~ machines are run at the T(5s) resonance, these methods could be 
used to investigate y in the B s . In the B s system, it has been estimated that y may be 
particularly large (5-15%) || and indeed similar methods have been previously discussed 
m@. 

In the cases of B s and D°, it is thought that CP violation in the mixing is small while 
in B°, mixing is known to produce large CP violating effects [JFL M. In our methods, CP 
violation in mixing will impact the signal, generally reducing it, so we will derive the 
formalism in a framework which includes CP violation. 
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The paper will proceed as follows. In Section 2 we will discuss the formalism which 
applies in the case where CP violation is negligible. In Section 3 we will discuss the 
application to D°D° production at a charm factory assuming there is no CP violation in 
D° oscillation or decay. In Section 4 we will generalize the formalism to the case where 
CP violation is present and consider the application to B and D mesons. In Section 5 
we will give our conclusions. In the appendix we give the time integrated correlated 
decay rates for H°H° decaying to various combinations of final states where indirect CP 
violation is present. 



2 Formalism if CP is Conserved 

The essential point of our method is best illustrated in cases where CP violation may 
be neglected. In such cases, p = q, so mass eigenstates also become eigenstates of CP, 
which we denote: 

\H ± ) = [\H°) ± \B°)] . (3) 

The crucial point here is that these CP eigenstates \H±) do not evolve with time. We 
can take advantage of this fact at threshold e + e~ machines, such as BaBar or CLEO-c. 
H°H° production at these machines is through resonances leading to HH pairs being in 
a quantum mechanically coherent state. Thus, if the production resonance has angular 
momentum L, the quantum mechanical state at the time of HH° production is 

V L = \H°H°) L = -L {iH^hWih)) + (-l) L \H°(k 2 )H°(h))} . (4) 

where k\ and k 2 are the momenta of the mesons. Rewriting this in terms of the CP basis 
we arrive at 

*2n = {\H + (k 1 )H + (k 2 )} - \H4h)H4k 2 )}} 

*2n + i = {{H-WH+fa)) + \H + (h)H4k 2 ))} (5) 

Thus in the L = odd case, which would apply to the experimentally important ^(3770) 
and T(4S) resonances, the CP eigenstates of the H mesons are anti-correlated while 
if L = even the eigenstates are correlatedF]. In either case the correlation between the 
eigenstates is independent of when they decay. In this way, if meson H(ki) decays to 
the final state which is also a CP eigenstate, then the CP eigenvalue of the meson H(k 2 ) 
can be determined. 



^^Whilc L = even resonances are not directly produced in e + e~ collisions, quantum-mechanically 
symmetric states can be produced in the decays, such as ^(4140) — > DD7. In the following, L = even 
case can also refer to this situation. 
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Using this eigenstate correlation as a tool to investigate CP violation has been sug- 
gested by §0 In this 

paper we will take advantage of such correlations for the experi- 
mental investigation of lifetime differences. The idea is fairly straightforward: we look at 
decays of the form ipi — * {H —* S a )(H — > Xlv) where S a is a CP eigenstate of eigenvalue 
a = ±1 and ipL generically means any resonance of angular momentum L that decays 
to HH° (see also a footnote after the Eqn. fl5|). 

In this scenario, the CP quantum numbers of the Hik^) is thus determined. The 
semi-leptonic width of this meson should be independent of the CP quantum number 
since it is flavor specific. It follows that the semi-leptonic branching ratio of iJ(A; 2 ) will 
be inversely proportional to the total width of that meson. Since we know whether Hik?) 
is a H + or and if_ from the decay of H{k\), we can easily determine y in terms of the 
semileptonic branching ratios of H±. 

This can be expressed simply by introducing the ratio 

r l = — •> (H — » S a )(H — » Xl^v)] 
CT T^ L ^(H^S a )(H^X)}Br(H°^Xlv)> [) 

where X in H —>■ X stands for an inclusive set of all final states. A deviation from 
R% = 1 implies a lifetime difference. From this experimentally obtained quantity, we 
extract y by 



3 Charmed Mesons if CP is Conserved 

Let us consider now the production of D°D° mesons at an electron positron collider. 
In the Standard Model, CP violation is expected to be small in D° hence the above 
formalism should apply directly to this case. For instance, the new tau-charm factory, 
under construction at CESR, will allow for simple and effective CP tagging in the case 
of -0(3770) — > DD because there are numerous candidates for S a in D° decay which have 
branching ratios in the few percent range, for instance Ks7t° (1.05%); K$uj (1.05%); 
K s rj' (0.85%); n+n- (0.15%). In additon, the modes K*°tt° and K*°p° may be used 
provided the K*° itself decays to a CP eigenstate, Ks,l^° and one can separate the main 
amplitude from cross channel processes. Similar comments apply to analogous states 
containing higher neutral kaonic resonances. 

These modes are thus candidates for S a in Eqn. @. If we write Eqn. (^) in terms of 
the semi-leptonic branching ratio of D±, B±, then equation Eqn. ([!]) becomes: 



\BUD) BUD)) y ' w 

2 For other measurements that involve CP correlations to study CP violation in Z)-mesons see [ [To[ . 
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In either case the statistical uncertainty in y is given by 



Ay = ( 2N A e B e A CT B a 



(9) 



where No is the initial number of ?/>'s, A a and A are the acceptances for the CP eigenstate 
modes and the semileptonic modes respectively while B a and B l are the branching ratios 
for those modes. In general, of course, we can combine the statistics for a number of 
modes so, as an example, if we assume that B e = 12%, B a = 2%, with A e A a = 0.1 then 
N = 10 8 gives Ay = 0.5%. 

At present, the information about the D° — D° mixing parameters x and y comes 
from the time-dependent analyses that can roughly be divided into two categories. First, 
more traditional analyses study time dependence of D — > / decays, where / is the final 
state that can be used to tag the flavor of the decayed meson. The most popular is the 
non-leptonic doubly Cabibbo suppressed (DCS) decay D — > K + tt~. Time-dependent 
studies allow one to separate the doubly Cabibbo suppressed decay from the mixing 
contribution, 



T[D°(t) -> K + ti- 



-rt I a 1 2 

\A K -TT+\ 



X 



R 2 



R + VRR m {y' cos - x' sin 0) Ft + (y 2 + x 2 ) (Ft) 2 



(10) 



where R is the ratio of Cabibbo favored (CF) and doubly Cabibbo suppressed decay 
rates and R m = \p/q\ while = arg(p/q). Since x and y are small, the best constraint 
comes from the linear terms in t that are also linear in x and y. Using this method, direct 
extraction of x and y is not possible from Eq. (|10J) due to unknown relative strong phase 
5 of DCS and CF amplitudes (see |TI| for extensive discussion), as x' = x cos 5 + ysin6, 
y cos 5 — x sin 5. This phase, however, can be measured independently [fT2 



The 



y __ 

corresponding formula can also be written for D° decay with x' — > —x' and R n 
0- 

Another method to measure D° mixing is to compare the lifetimes extracted from 
the analysis of D decays into the CP-even and CP-odd final states. This study is also 
sensitive to a linear function of y, via 



t(D K+K- 



1 = y cos + x sin < 



R 2 



(11) 



Time-integrated studies of the semileptonic transitions are sensitive to the quadratic 
form x 2 + y 2 and at the moment are not competitive with the analyses discussed above. 



4 Formalism if CP is Violated in H°H° Oscillations 



In the case where CP violation is present in the H°H° mixing, it is necessary to consider 
general time dependent entangled states of the H°H° pair. Following the notation of 



4 



we will denote the wave function \H(t)) at a given moment in time, t, by a two element 
vector: 

\H(t)) = a\H°) + a\H°)= (12) 
CPT conservation forces the general mass matrix in the following form 




M = M + iT/2 =\ (13) 



where M and f are Hermitian while A, p and q are in general complex numbers. The 
effects of CP violation in the system are usually parameterized in terms of the ratio: 

V - = \±± = Rm j* (14) 
q 1 — e 

In the limit of CP conservation in mixing matrix, R m = 1. Even if CP is violated, in 
the case of heavy neutral mesons, it is expected that R m ~ 1. The phase 0, of course, 
depends on the convention one uses for weak phases that can be traded off against the 
weak phase in the decay in the usual way. In our discussion it will be useful to assume 
that we are using a convention where we have absorbed any weak phase from the decay 
into the the mixing. 

For an isolated H meson, the wave function at time t is related to the the wave 
function at t — by: 

\H(t)) = U t \H(0)), (15) 
where the time evolution operator U t satisfies the equation 

i^=MU u U = l. (16) 
This Schrodinger-like equation can be solved to yield the familiar result 

u = ( 9+{t) Wtid-M] (i 7 ) 
1 \(q/p)9-(t) 9 + (t) )' 1 ) 

Here, the time dependence of D° and D° is driven by 

g+(t) = (cosh yr/2 cos xt/2 — i sinhyr/2 sin xt/2) e~ Mr / 2 

g-(t) = (-sinh yr/2 cos xr/2 + i cosh yr/2 sin rrr/2) e~ Mr/2 (18) 

with t — Ft and /i — 1 + 2im/Y. 
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Let us now consider the time integrated decay rate for a single if to a final state /. 
If a and a are the amplitudes for H° and H° to decay to / respectively and l^o) is the 
initial wave function for the meson, then the time integrated decay rate is 

2I70o) = (Q + P)tr[p fPo } + (Q- P)tr[v^p fP o] - 2Re [(yQ - ixP)tr[p f vp )] , (19) 

where 

and 

P = l/(l + a; 2 ) Q=l/(l-y 2 ). (21) 

Of particular interest is the where / is a CP eigenstate with CP=a = ±1. If 

we assume R m = 1, as would be the case for E>d — ► ipK s in the Standard Model, then 



l(r 1 (H°) + r l (H°))= 1 -^^r (22) 



y 

l(r f (H°)-r,m) = ^r«, (23) 

where T is the decay rate for H to /. 

This can be easily generalized to the case of the entangled initial state which presents 
itself in the creation of H°H° pairs from a ^ resonance. As was shown in Eqns. (^HJ), 
the states of interest can be decomposed into the coherent sum of products of flavor (or 
CP) eigenstates. Using Eqn. (Q) we can write the time integrated correlated decay rate 
for V L - (H - f x ){H -/ 2 ) is: 

r /l/2 (^) = r fl (H°)T f2 (H°) + r h (H°)T h (H ) 
+(-i) i [r /l (p + -)r /2 (p- + ) + r /l (p_ + )r /2 (p + -)] (24) 

where p^ are the matrices: 



P++ = \ n „ P 





p+- = L J p-+ = L J ( 25 ) 



Let us return to the calculation of y through the determination of R% defined in Eqn. @. 
In the case L = 1 the numerator is F s " Xl (^fi) and so Eqn. (pl|) implies: 



^2 + * 2 (i^ 2 + l)+y 2 (P± 2 -i; 



1 L) 2(l + x 2 )(l-y 2 ) 1 1 
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where we assume there is no further CP violation in the decay amplitude, and the ± 
signs are for the positively and negatively charged leptons respectively. 

The denominator of Eqn. (|^) is given by Eqn. (p4|) where f\ is 5*0- and all possible 
values of fi are summed over. In this case the L dependent terms vanish and the rest 
simplifies to: 



h 



(27) 



where 1 is the identity matrix. This is proportional to the average between the time 
integrated decay rate of H° and H° to the final state f\. It is easiest to see that in the 
CP-eigenstate basis spanned by the states of Eq. (|5]). Indeed, in the particular example 
of L = odd the time-integrated decay rate is 



T flX (^ 1 ) 



dtidt 2 - 



Tr 

+ Tr 
— Tr 



Ut2 Px U t2 P-h 
Ut2 PX U t2 P— 

U ] t2 p x U t2 P+- 



Tr 

Tr 
Tr 



Uti Pn u n P- 
U\x p n U t i Ph 

u\i Pn Un p- 



(28) 



- Tr [ U} 2 p x U t2 p-+] Tr [ p fl U tl p + _ 



where px and pfi are the matrices for H — >• X and H —> fi decay amplitudes respec- 
tively (see Eqn (pO|)). It is easy to see that in the mass eigenstate basis ll} 2 Px U t2 \mass — 
diag(Tie~ rit2 , r 2 e _r2 * 2 ). In principle, this needs to be translated back to the CP eigen- 
state basis. However, integration with respect to t 2 yields a unit matrix, which is in- 
variant under the change of basis. This simplifies the Eq. ( p8j) considerably, which, after 
taking the corresponding traces transforms into 



r^v^) = ~ J dt Tr [pfiUtUj 

which, in the limit R m = 1, becomes for the semileptonic final state (a complete expres- 
sion is available in the Appendix) 



T XX '(V L ) oc 



1 + (-lr cos. 



(i-y 2 ). 

In the limit R m = 1 the ratio of Eqs. (PB[) and (B^) becomes: 



1 + (— l) L <jy cos . 
In which case the generalization of Eqn. (0) is 



ycos. 



(29) 



(30) 



(31) 
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So we can regard the measurement of R% as leading to a determination of y cos 0. A 
similar result holds for the non-leptonic final state (such as D — > Ktt, with corrections 
proportional to R) . In the case where R m ^ 1 the corresponding expression depends on 
x as well as y. For instance, for L = 1 

^ L PQ(1 + P m x 2 cosh a m + P m y 2 sinh a m ) 

(Q cosh 2 cx m — P sinh 2 a m — xP sinh a m sin — yQ cosh a m cos 0) 

where a m = log(P m ) = log \/l + A m « A m /2 fli~3"| . Expanding this to first order in a m 
we obtain: 



pL 1 (x 2 + y 2 )(l -ycos0) +x(l -y 2 )sin0 , 2 , , v 

1 — cry cos (1 — y cos 0)^(1 + x l ) 

Thus, if we define y by 

ycos0 = cr a RL (34) 



If we expand y to first order in a m we obtain: 



y = y 



(x 2 + y 2 )(l — y cos0) + (1 — y 2 )xsi 



sm i 



[1 + X 2 ) COS( 



(35) 



Clearly then, Eqn. ( |3~1~D gives y only if a m is known to be small. The actual value of a m 
can be experimentally obtained from the semileptonic decay asymmetry ||. 

In our discussion we will now assume that R m ~ 1 and so the ratio R^ gives us 
ycos0 through Eqn. ([H]). The error in determining y is thus given by the generalization 
of Eqn. (BD 



Ay cos 0= (2N A £ B £ A a B a ) ' (36) 

In the case of Z)°, the systematics for Ay cos is the same as the systematics for Ay in 
the CP conserving case discussed above. 

In the case of B°, which is equal to 2(3 in the Standard Model has been measured 
at the BaBar and BELLE experiments |7], |) The average of these two results is currently 
sin 2(3 = 0.78±0.08 thus cos 2/3 ~ 0.6. Thus, if we take N Q = 10 8 and use only ipKs decay 
mode with ip — ► and assume that A a A\ ~ 1/4 then Ay cos 2(3 = 0.06 corresponding 
to Ay = 0.1. Clearly bringing in additional 5*0- modes will improve determination of 
Ay. We can also improve the statistics by using flavor specific decays of the B° other 
than pure leptonic decays. The BaBaR and Belle experiments have made considerable 
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progress in their ability to accomplish this and obtain an effective value of AiB\ pa 0.7. 
Using this result, the above gives Ay pa 0.026. 

To produce correlated B s pairs one needs to run an electron-positron machine at the 
T(5s) resonance. This state can decay into B S B S , B*B S and B*B* where the B* decays 
to B s j. As discussed in || if there are or 2 photons in the final state (i.e. the decay 
was to B S B S or B*B*) then the B S B S is in an L = odd state while if there is one photon 
in the final state (i.e. B*B S ) then the final B S B S state is L = even. 

The branching ratio to S a states in the case of B s is in principle much larger than 
in the case of Bq. For instance, the branching ratio for B s — > DfD~ should be similar 
to the measured branching ratio for B° — > D~D+ which is about 0.8%. likewise one 
can also estimate the branching ratio of B s — > J/^rf'^ at about 0.15% in addition 
analogous states such as D*D* etc should have branching ratios on the order of 0.1% at 
least. The acceptance for such states may be lower than for ipK s so we will assume that 
A„Ai ~ 0.1 with a branching ratio to CP states of about 0.8%. Using these assumptions, 
if one had a high luminosity T(5s) machine that was able to produce 10 8 B s pairs then 
Ay cos = 0.7% which would be the same as Ay if the Standard Model expectation 
that = was correct. In |J it was shown that using a method of generalizing the 
identification of CP eigenstate decays of B s to include all states with a final quark content 
of cess at the expense of the efficiency of such tagging may allow the determination of 
Ay to a precision of about 0.28% under the same assumption although in this case the 
precision also depends on the value of y. 



5 Conclusions 

In summary, we discussed the possibility of time-independent measurements of lifetime 
differences in D and B systems. It is important to reiterate that time-dependent mea- 
surements are quite difficult at the symmetric e + e~ threshold machines due to the fact 
that the pair-produced heavy mesons are almost at rest [14]]. The techniques described 



above will provide a time integrated quantity that is separately sensitive to the lifetime 
difference y. 

This will be particularly useful in the case of D° where a charm factory running at 
the ■0(3770) resonance can yield the measurement with precision of Ay cos <fi ~ 0.5% 
which is in the range of some standard model predictions. At a T(5S) B factory with a 
luminosity sufficient to produce 10 8 B s pairs, a precision of Ay cos <fi pa 0.7% should be 
achievable which is much smaller that the standard model prediction of 5-15%. Thus, 
even if only ~ 10 6 B s pairs are produced, precision on the order of the Standard Model 
prediction can be obtained. In the case of the T(4S) B factory with 10 s BB pairs, 
Ay cos 2(3 pa 3% which does not probe to the level of the Standard Model estimate in 
this case. Yet, a new high luminosity B factory (such as discussed SuperBaBar) will be 
able to measure y in B d system. 
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Appendix: Correlated Decays with CP Violation 

In this appendix we provide the expressions for the time integrated decay of a correlated 
H° -TP state to various pairs of final states using the formalism discussed in the text. 
The final states we consider are: 

1. S±, a CP eigenstate such as in D° — > K s n° or B d ( s ) — > J/i/jK s (4>). 

2. L ± ) a flavor specific semi-leptonic decay to a final state containing £ ± . 

3. G, a hadronic final state such that both H° and H° can decay to it. For example, 
in charmed mesons, D° — > G + is Cabibbo Favored (CF) and D° — > G~ is doubly 
Cabibbo suppressed, as in D° — > K~n + . This implies that the ratio of DCS to 
CA decay rates R is small and the results can be expanded in terms of this ratio. 
Alternatively, the ratio of amplitudes can be of order one in B decay, as in the 
example of B s — > DfK~, so all powers of R must be kept. 

4. X is an inclusive set of all final states. 

It is now easy to construct all possible combinations of the above final states. For the 
case of antisymmetric initial state (L =odd), we have for T^j 2 



odd 



2 + {l + R^)x 



r (S a )T (L+) 
2(l-|/ 2 )(l + x 2 ) 



(37) 



r s+s. 



8Ri + (l + Ri)(x 2 + y 2 ) 

+2R 2 m ((l + 2 cos 2 (f))x 2 + (1 + 2 sin 2 (f))y 2 ) 
r (S + )T (S-) 



x 



2(1 -y 2 )(l + x 2 ) 



x 



(38) 



odd 



(x 2 + y 2 )(cosh 2 a m — cos 2 



[l-y 2 )(l + x 2 ) 



(39) 
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r^ L± = {R m \x 2 + y 2 ) 



2(1 - y 2 )(l + x 2 ) 



odd 



2 + x 2 - y 2 



2(1 - y 2 )(l + x 2 ) 



odd 



1 + y 2 sinh 2 a m + x 2 cosh 2 a m 



—<j(y{l + x 2 ) cosha m cos0 + x(l — y 2 ) sinha m sin 
2T D T (S a ) 



X 



(1 -y 2 )(l + x 2 ) 



x 



r L CT X _ 
1 odd — 



1 + x 2 cosh 2 a m + y 2 sinh 2 a m H — (x 2 + y 2 ) sinh 2a r , 

2 



X 



X 



2r D r (L+) 

(l_ J/ 2)( 1+a .2) 



odd 



+ y 2 +(l + R 2 )(2 + x 2 - y 2 ) R 2 m + R 2 (x 2 + y 2 ) R 
+2Rx (l - y 2 ) sin(5 + 4>)R m (l - R 2 m ) 
-2Ry (l + x 2 ) cos(5 + (f))R m (l + R 

T(G + )T D 
A{l-y 2 ){l+x 2 )R 2 m 



-pGScr 
1 Odd 



R 2 m {\ + R 2 ) + (R 2 + R 2 J(1 + R 2 Jx 2 + (R 2 - R 2 J(1 - R 2 Jy 2 
+4ri? m (y 2 (cos0sin0sin<5 — sin 2 cos 5) 
+x 2 (cos sin sin 5 + cos 2 cos 5) + cos 
2T(G + )r(S a ) 



x 



x 



(1 -y 2 )(l + x 2 
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odd 



2R 2 + (Rl 2 + R 2 )x 2 + (Rl 2 - R 2 )y' 



r(G)r(L+) 

2(l-y 2 )(l+x 2 



(46) 



odd 



2 + (R'R% + l)x 2 + (R 2 R±? - l)y 



r(G)r(L- 



2(1 -y 2 )(l +x 2 ) 



(47) 



<?±G=F 
odd 



(l + i? 4 ) i& (-2 - x 2 + y 2 ) - R 2 (x 2 + y 2 



R 2 Ri (x 2 + 



+ 2R 2 R 2 



(2 + x 2 - y 2 ) cos 25 + (x 2 + y 2 ) cos : 

r 2 (G) 



x 



2(y 2 -l)(l + x 2 )R 2 m 



(48) 



,G±G± 
odd 



- 2 + y 2 ){R 4 R^ - 2R 2 cos 2(5 ± 0) + i?^ 2 ) 



x 



x 



r 2 (<^ 



2(1 -y 2 )(l +x 2 ) 



(49) 



In the case of charmed mesons, R <C 1. Neglecting possible CP-violating effects 
and taking the ratio of Eqs. (fl~9|) and (|48|) simultaneously expanding numerator and 
denominator in R, x, and y we reproduce the well-known result that DCS/CF interference 
cancels out in the ratio for L =odd |15| and gives the result, identical to the semileptonic 
final state, (x 2 + y 2 )/2. 

The results for L =even are more cumbersome, so we present only a few of r{^: 



(x 2 + y 2 ) (3+(x 2 -y 2 )+x 2 y 2 ) 
+Ri 2 (2 + (1 + x 4 )x 2 - (l - 4x 2 - x 4 ) y 2 + (l - x 2 ) y 4 ) 



-AR^a ( (l + £ 2 ) 2 y cos <p =F (l - y 2 ) 2 x sin . 



r (g g )r (L+) 

2(l-y 2 ) 2 (l + x 2 ) 2 



(50) 



x 2 + y 2 )(x 2 + (x 2 - l)y A + 3)(cosh2a m - cos 20) 



(l-y 2 ) 2 (l + x 2 y 



(51) 



,L±L± 



^ 2 (^ 2 + y 2 )(^ 2 + (x 2 -l)l/ 2 + 3) 



r 2 (^ H 



2(l-?/ 2 ) 2 (l + x 2 ) 2 



(52) 
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x 4 + x 2 - (x 2 - l)y 4 + (x 4 + Ax 2 - l)y 2 + 2) 



2(1 -?/ 2 ) 2 (l + x 



,,2^2 



(53) 



Taking the ratios of the decay rates presented above, one can easily generalize the 
results of [IS| to the case of CP-nonconservation. 
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